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The first practical method to evolve many-body nuclear forces to softened form using the Similarity 
Renormalization Group (SRG) in a harmonic oscillator basis is demonstrated. When applied to *He 
calculations, the two- and three-body oscillator matrix elements yield rapid convergence of the 
ground-state energy with a small net contribution of the induced four-body force. 

PACS numbers: 21.30.-x,05.10.Cc,13.75.Cs 



A major goal of nuclear structure theory is to make 
quantitative calculations of low-energy nuclear observ- 
ables starting from microscopic internucleon forces. Chi- 
ral effective field theory (xEFT) provides a systematic 
construction of these forces, including a hierarchy of 
many-body forces of decreasing strength Renormal- 
ization group (RG) methods can be used to soften the 
short-range repulsion and short-range tensor components 
of the initial chiral interactions so that convergence of nu- 
clear structure calculations is greatly accelerated [1, [3] . 
The difficulty is that these transformations (or any other 
softening transformations) change the short-range many- 
body forces. To account for these changes, we present 
in this letter the first consistent evolution of three-body 
forces by using the Similarity Renormalization Group 
(SRG) which offers a technically simpler 

approach to evolving many-body forces than other RG 
formulations. Our results show that both the many-body 
hierarchy of ^EFT and the improved convergence prop- 
erties are preserved. 

The SRG is a series of unitary transformations of the 
free-space Hamiltonian {H = H\=oc), 

Hx = UxHx=ooUi , (1) 

labeled by a momentum parameter A that runs from oo 
toward zero, which keeps track of the sequence of Hamil- 
tonians (s = has been used elsewhere 0, Si)- These 
transformations are implemented as a flow equation in A 
(in units where Y? jM = 1), 

f^ = 4r,i/.],i/.], (2) 

whose form guarantees that the are unitarily equiv- 
alent [1,0]. 

The appearance of the nuclcon kinetic energy T in 
Eq. ^ leads to high- and low-momentum parts of ii\ 
being decoupled, which means softer and more conver- 
gent potentials [9| . This is evident in a partial- wave mo- 
mentum basis, where matrix elements (lii\H\\k') connect- 
ing states with (kinetic) energies differing by more than 
A^ are suppressed by e~^^^~^ ^)^/a'* factors and there- 
fore the states decouple as A decreases. (Decoupling also 
results from replacing T in Eq. ^ with other genera- 
tors [i, 0, 0, [ill-) The optimal range for A is not yet 



established and also depends on the system, but experi- 
ence with SRG and other low-momentum potentials sug- 
gest that running to about A = 2.0 fm^^ is a good com- 
promise between improved convergence from decoupling 
and the growth of induced many-body interactions [(!]. 
(Also, differences between using T and the diagonal of 
H\ in Eg. d^, which can be very important in some sit- 
uations [iO|, are negligible in this A range.) 

To see how the two-, three-, and higher-body potentials 
are identified, it is useful to decompose H\ in second- 
quantized form. Schematically (suppressing indices and 
sums), 

Hx = {T)a)a -f {V^^^)a)a)aa + {V^^'^)a^ a) aaa , 

(3) 

where a\ a are creation and destruction operators with 
respect to the vacuum in some (coupled) single-particle 

basis. This defines (T), {V^^^), (V^^^^), ...as the one- 
body, two-body, three-body, . . . matrix elements at each 
A. Upon evaluating the commutators in Eq. ([2]) using 
Hx from Eq. ([3]), we see that even if initially there are 
only two-body potentials, higher-body potentials are gen- 
erated with each step in A. Thus, when applied in an 
A-body subspace, the SRG will "induce" A-body forces. 
But we also see that (T) is fixed, (V^^^^) is determined 
only in the A = 2 subspace with no dependence on (V^^ ) , 
(V^^) is determined in A = 3 given {V^^), and so on. 

Since only the Hamiltonian enters the SRG evolution 
equations, there are no difficulties from having to solve 
T matrices in all channels for different A-body systems. 
However, in a momentum basis the presence of spectator 
nucleons requires solving separate equations for each set 
of {V^'') matrix elements. In Refs. a diagram- 

matic approach is introduced to handle this decomposi- 
tion. But while it is natural to solve Eq. ([2]) in momentum 
representation, it is an operator equation so we can use 
any convenient basis. Here we evolve in a discrete basis, 
where spectators are handled without a decomposition 
and induced many-body forces can be directly identified. 
Having chosen such a basis, we obtain coupled first-order 
differential equations for the matrix elements of the flow- 
ing Hamiltonian Hx, where the right side of Eq. ^ is 
evaluated using simple matrix multiplications. 
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Our calculations are performed in the Jacobi coordi- 
nate harmonic oscillator (HO) basis of the No-Core Shell 
Model (NCSM) [H]. This is a translationally invariant, 
anti-symmetric basis for each A, with a complete set of 
states up to a maximum excitation of Nj^i^^hfl above the 
minimum energy configuration, where f2 is the harmonic 
oscillator parameter. The procedures used here build di- 
rectly on Ref. J 3.] , which presents a one-dimensional im- 
plementation of our approach along with a general anal- 
ysis of the evolving many-body hierarchy. 

We start by evolving H\ in the A = 2 subsystem, which 
completely fixes the two-body matrix elements (V^ ). 
Next, by evolving H\ in the A = S subsystem we deter- 
mine the combined two-plus-three-body matrix elements. 
We can isolate the three-body matrix elements by sub- 

(2) 

tracting the evolved {V^ } elements in the A = 3 ba- 
sis [l^. Having obtained the separate NN and NNN 
matrix elements, we can apply them unchanged to any 
nucleus. We are also free to include any initial three- 
nucleon force in the initial Hamiltonian without chang- 
ing the procedure. If applied to A > 4, four-body (and 
higher) forces will not be included and so the transforma- 
tions will be only approximately unitary. The questions 
to be addressed are whether the decreasing hierarchy of 
many-body forces is maintained and whether the induced 
four-body contribution is unnaturally large. We summa- 
rize in Table H] the different calculations to be made for 
■^H and ''He to confront these questions. 

The initial (A = oo) NN potential used here is the 
500 MeV N^LO interaction from Ref. The initial 

NNN potential is the N^LO interaction in the local 
form of Ref. [l3| with constants fit to the average of tri- 
ton and ^He binding energies and to triton beta decay 
according to Ref. We expect similar results from 

other initial interactions because the SRG drives them 
toward near universal form; a survey will be given in 
Ref. [l^. NCSM calculations with these initial interac- 
tions and the parameter set in Table I of Ref. ^Ts'l yield 
energies of -8.473(4) MeV for and -28.50(2) MeV for 
'^He compared with -8.482 MeV and -28.296 MeV from 
experiment, respectively. So there is a 20 keV uncertainty 
in the calculation of ^He from incomplete convergence 
and a 200 keV discrepancy with experiment. The latter 
is consistent with the omission of three- and four-body 
chiral interactions at N'^LO. These provide a scale for 
assessing whether induced four-body contributions are 
important compared to other uncertainties. 

In Fig. [H the ground-state energy of the triton is plot- 
ted as a function of the flow parameter A. Evolution is 
from A = oo, which is the initial (or "bare") interaction, 
toward A = 0. We use TV^ax = 36 and Hfl = 28 MeV, for 
which all energies are converged to better than lOkeV. 
We first consider an NN interaction with no initial NNN 
("NN-only"). If H\ is evolved only in an A = 2 sys- 
tem, higher-body induced pieces are lost. The resulting 
energy calculations will only be approximately unitary 
for A > 2 and the ground-state energy will vary with A 
(squares). Keeping the induced NNN yields a flat line 
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FIG. 1: (Color online) Ground-state energy of as a func- 
tion of the SRG evolution parameter, A. See Table |T] for the 
nomenclature of the curves. 




FIG. 2: (Color online) Ground-state energy of *He as a func- 
tion of the SRG evolution parameter, A. See Table |T] for the 
nomenclature of the curves. 



(circles), which implies an exactly unitary transforma- 
tion; the line is equally flat if an initial NNN is included 
(diamonds). Note that the net induced three-body is 
comparable to the initial NNN contribution and thus is 
of natural size. 

In Fig. [21 we examine the SRG evolution in A for ^He 



with hn = 36 MeV. The (V^^^) and {V^"') matrix ele- 
ments were evolved in ^ = 2 and A — S with iVmax = 28 
and then truncated to A^max = 18 at each A to diagonal- 
ize ''He. The NN-only curve has a similar shape as for 
the triton. In fact, this pattern of variation has been ob- 
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TABLE I; Definitions of the various calculations. 

NN-only No initial NNN interaction and do not keep NNN-induced interaction. 

NN + NNN-induced No initial NNN interaction but keep the SRG-induced NNN interaction. 

NN + NNN Include an initial NNN interaction and keep the SRG-induced NNN interaction. 



served in all SRG calculations of light nuclei [3]. When 
the induced NNN is included, the evolution is close to 
unitary and the pattern only depends slightly on an ini- 
tial NNN interaction. In both cases the dotted line rep- 
resents the converged value for the initial Hamiltonian. 
At large A, the discrepancy is due to a lack of conver- 
gence at A'^inax = 18, but at A < 3fm~^ SRG decouphng 
takes over and the discrepancy is due to short-range in- 
duced four-body forces, which therefore contribute about 
50keV net at A = 2fm~^. This is small compared to the 
rough estimate in Ref. [2^ that the contribution from 
the long-ranged part of the N^LO four-nucleon force to 
^He binding is of order a few hundred keV. If needed, 
we could evolve 4-body matrix elements in ^ = 4 and 
will do so when nuclear structure codes can accomodate 
them. 
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FIG. 3: (Color online) Ground-state energy of H as a func- 
tion of the basis size A^max for an N^LO NN interaction [l^ 
with and without an initial NNN interaction P, [l^. Un- 
evolved ("bare") and Lee-Suzuki (L-S) results with hfl = 
28MeV are compared with SRG at hfl = 20MeV evolved 
to A = 2.0 fm-^ 

In Fig. [3l we show the triton ground-state energy as 
a function of the oscillator basis size, A^max, for various 
calculations. The lower (upper) curves are with (with- 
out) an initial three-body force (see Table [!}. The con- 
vergence of the bare interaction is compared with the 
SRG evolved to A = 2.0 fm"^ The oscillator parame- 
ter hn in each case was chosen roughly to optimize the 
convergence of each Hamiltonian. (As A decreases, so 
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FIG. 4: (Color online) Ground-state energy of *He as a func- 
tion of the basis size Amax for an N^LO NN interaction 
with an initial NNN interaction P, [ll]. Unevolved (bare) re- 
sults are compared with Lee-Suzuki (L-S) and SRG evolved 
to A = 2.0 fm"^ at HQ = 28 and 36MeV. 



does the optimal h^l.) We also compare to a Lee-Suzuki 
(L-S) effective interaction, which has been used in the 
NCSM to greatly improve convergence [U, These 
effective interactions result from unitary transformations 
within the model space of a given nucleus, in contrast to 
the free-space transformation of the SRG, which yields 
nucleus-independent matrix elements. 

The SRG calculations are variational and converge 
smoothly and rapidly from above with or without an 
initial three-body force. The dramatic improvement in 
convergence rate compared to the initial interaction is 
seen even though the xEFT interaction is relatively soft. 
Thus, once evolved, a much smaller A'max basis is ade- 
quate for a desired accuracy and extrapolating in A'max 
is also feasible. 

Figure |3] illustrates for ^He the same rapid conver- 
gence with A'max of an SRG-evolved interaction. How- 
ever, in this case the asymptotic value of the energy dif- 
fers slightly because of the omitted induced four-body 
contribution. (The SRG-evolved asymptotic values for 
different h^l differ by only lOkeV, so the gap between 
the converged bare/L-S results and the SRG results is 
dominated by the induced NNNN rather than incomplete 
con verg ence) . Convergence is even faster for lower A val- 
ues [lOj, ensuring a useful range for the analysis of few- 
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FIG. 5: (Color online) Binding energy of the alpha particle 
vs. the binding energy of the triton. The Tjon line from phe- 
nomenological NN potentials (dotted) is compared with the 
trajectory of SRG energies when only the NN interaction is 
kept (circles). When the initial and induced NNN interac- 
tions are included, the trajectory lies close to experiment for 
A > 1.7 fm~^ (see inset). 



analogous to the well-known Tjon line (dotted), which 
is the approximate locus of points for phenomenological 
potentials fit to NN data but not including NNN [23] . In 
contrast, the short trajectory of the SRG with the NN 
+ NNN interaction (shown for A > 1.8 fm~^) highlights 
the small variations from the omitted four-nucleon force. 
Note that a trajectory plotted for NN-f NNN-induced cal- 
culations would be a similarly small line at the N'^LO 
NN-only point. 

In summary, we have demonstrated a practical method 
to use the SRG to evolve NNN (and higher many-body) 
forces in a harmonic oscillator basis. Calculations of 
A < A nuclei including NNN show the same favorable 
convergence properties observed elsewhere for NN-only, 
with a net induced four-body contribution in A — A that 
is smaller than the truncation errors of the chiral inter- 
action. The soft SRG interactions are an alternative to 
the use of Lee-Suzuki effective interactions in NCSM and 
the HO matrix elements can also be used (after conver- 
sion to a Slater-determinant HO basis as needed) for cou- 
pled cluster and many-body perturbation theory calcu- 
lations. A more complete analysis of convergence and 
dependencies for the energy and other observables for 
few-body systems, as well as results for other interac- 
tions and choices of generator in Eq. ([2]), will be given in 
a forthcoming publication p^ . 



body systems. However, because of the strong density 
dependence of four-nucleon forces, it will be important 
to monitor the size of the induced four-body contribu- 
tions for heavier nuclei and nuclear matter. 

The impact of evolving the full three-body force is 
neatly illustrated in Fig. [5l where the binding energy of 
"'He is plotted against the binding energy of ^H. The ex- 
perimental values of these quantities, which are known 
to a small fraction of a keV, define only a point in this 
plane (at the center of the X, see inset) . The SRG NN- 
only results trace out a trajectory in the plane that is 
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